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Most fully developed galaxies have a vivid spiral structure, but the formation and 
evolution of the spiral structure are still an enigma in astrophysics. In this paper, 
according to the standard Newtonian gravitational theory and some observational 
facts, we derive an idealized model for spiral galaxy, and give a possible explanation 
to the spiral structure. We solve some analytic solutions to a spiral galaxy, and obtain 
manifest relations between mass density and speed. Prom the solution we get some 
interesting results: (I) The spiral pattern is a stationary or static structure of density 
wave, and the barred galaxy globally rotate around an axis at tiny angular speed. 
(II) All stars in the disc of a barred spiral galaxy move in almost circular orbits. (Ill) 
In the spiral arms, the speed of stars takes minimum and the stellar density takes 
maximum. (IV) The mass-energy density of the dark halo is compensatory for that 
of the disc, namely, it takes minimum in the spiral arms. This phenomenon might 
reflect the complicated stream lines of the dark halo. 

PACS numbers: 98.52.Nr, 98.62.Hr, 98.62.Dm, 98.10,+z, 98.35.Df 



I. INTRODUCTION 

Most fully developed galaxies have a vivid spiral structure, which has profound influences 
on the generation and evolution of stars. During the last 60 years many efforts have been 
made to reveal the nature of the spiral structure. Despite the progress made in understanding 
the amplification mechanisms of the spiral arms and the role of the spiral arms in the 
dynamics of galaxies, the formation and evolution of the spiral structure are still an enigma 



*Electronic address: yqgu@fudan.edu.cn 



2 



in astrophysics [1, 2, 3]. We even do not know if the spirals are a long-lived phenomenon, or 
are short-lived and regenerated many times during the galactic evolution. C. C. Lin and F. 
H. Shu suggested that, the spiral structure should be a quasi-stationary and stable density 
waves, which are excited by over-reflection at co-rotation due to some nonlinear effects[4, 5]. 
But some other opinions favor short-lived or recurrent spiral patterns which are developed in 
the galactic disks via swing-amplification in a rotating disk with shear [6], or by an external 
force[7]. The recent observational results seem to favor a long-lived pattern. In contrast 
the nearby galaxies with the distant ones observed through the Hubble Space Telescope and 
from the ground[8, 9], the results show that there is only mild evolution in the relationship 
between radial size and stellar mass for galactic disks within the redshift region < z < 1. 

In the noncentral potential, the orbit of a star is usually an unclosed and complicated 
spatial curve, which sensitively depends on the initial speed and direction[l]. So the stellar 
system should be more reasonably treated as fluid rather than mass points. Although the 
dynamics for the two models is essentially the same Newtonian mechanics, but the initial 
and boundary conditions are different. In the case of fluid, the stream lines are consistent 
field, which are a natural result of the generation of stars from nebulae. 

W. Dehnen and J. Binney fitted the mass distribution within the Milky Way to the 
observational data such as the rotation curve and Oort's constants[10], but they found 
the fitting mass distribution is ill determined. In [11, 12], the authors made numerical 
simulations for 2-dimensional stellar hydrodynamics of a flat galactic disc embedded in dark 
matter halo, and the model solves the Boltzmann equations up to second order moment [1]. A 
dynamical approach to explain the formation of both spirals and rings in barred galaxies was 
proposed in [13, 14]. It is based on the orbital motion driven by the unstable equilibrium 
points of a given rotating bar potential, and then the spirals, rings and pseudo-rings are 
related to the invariant manifolds associated to the periodic orbits around these equilibrium 
points. By adjusting dynamical parameters of the host galaxy, we get the spiral and ring 
like structures. 

If arms of a spiral galaxy were constructed by fixed material, the arms would become 
more and more tightly wound, since the matter nearer to the center of the galaxy rotates 
faster than the matter at the edge of the galaxy. The arms would become indistinguishable 
from the rest of the galaxy after only a few orbits. In order to avoid this winding problem, 
Lin and Shu proposed the density wave theory [4, 5]. They suggested that the arms are not 
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constructed by stationary materials, but instead made up of areas of greater density similar 
to a traffic jam on a highway. More specifically, for some dynamical reasons, the luminous 
materials such as stars intend to move slowly and spend more time inside the spiral arms 
and move rapidly outside the spiral arms. Then an inhomogeneous and quasi-stationary 
distribution of mass density wave is developed, which forms the spiral arms. In the sense 
of morphology, this theory is quite successful in some aspects. However, the dynamical 
explanation for the formation of density wave is still unclear. The calculations of this paper 
attempt to give a clear dynamical explanation for the density wave. 

According to the following observational facts, we find the problem can be simplified and 
described by an idealized fluid model, and then the nature of the spiral structure might be 
disclosed by dynamical approach. 

1. By observation, except for the center part, all stars in a spiral galaxy are mainly 
distributed in a thin disc, so the spiral structure should be a visible display of collective 
motion of these stars, and the collective motion can be approximately described by 
the 2-dimensional hydrodynamics. 

2. The collision among stars rarely occurs, so the fluid of star flow is zero-pressure and 
inviscid, and the stars move along geodesies. 

3. All stars are driven by the average background gravity of total mass-energy distribution 
in the galaxy, and the background mass density and gravity can be constrained by some 
empirical data such as the flat rotation curve. 

These are some basic assumptions of the following calculations. The results show that, the 
stationary solutions to such system indeed exist, and the distributions of mass density and 
speed do have stable spiral structure. So this idealized model and the solutions might shed 
some lights on the enigma of galactic structure. 

II. HYDRODYNAMICS FOR STARS IN A SPIRAL GALAXY 

As discussed below, a fully dynamical approach to the galactic structure in the context 
of general relativity strongly depends on the properties of the dark halo, so it is unrealistic 
at present due to the lack of such knowledge[15, 16]. Since the gravity in a galaxy is very 



i 



weak except for the region near the center, the Newtonian gravitational theory is accurate 
enough to simulate the spiral structure. 

Denote the total effective mass-energy density of a galaxy by p, the Newtonian gravita- 
tional potential by $, we have the dynamical equation for Newtonian gravity (see [17] or 
the appendix) 

d a d a § = -AirGp, (2.1) 

where d a d a = d 2 — V 2 is the d'Alembert operator. For the stars moving in the gravity, the 
dynamical equation of the flow is given by 

(d t + V ■ V)V = -V$. (2.2) 

Theoretically, the complete dynamical equation system should include the continuity equa- 
tion and the dynamics of the dark halo. But these equations strongly depend on the equation 
of state of dark matter and dark energy, which is not clearly discovered. However, as shown 
bellow, the lack of this knowledge can be compensated by other empirical data such as the 
flat rotation speed curve, and then some important information of galaxy such as the total 
mass density and the spiral structure, can be inversely derived. That is to say, we do the 
inverse treatments as done in [10, 13, 14]. 

For convenience we rewrite (2.1) and (2.2) in the spherical coordinate system (t,r,9,<j)). 
By 3-dimensional tensor operation, we get the following equations 

c~ 2 d 2 <5> - [d 2 r + -d r + \{d 2 + cot Ode + — U35)]$ + 47rGp = 0, (2.3) 
r r A sin 9 v 

(d t + V r d r + V e d e + V^)V r - rV e 2 - r sin 2 6V 2 + d r $ = 0, (2.4) 
(dt + V r d r + V e de + V^)V + \v e - sin 9 cos 6V 2 + ^ = 0, (2.5) 

(dt + V r d r + Vede + V^Vj + -V r V^ + 2 cot 9V B V <t> + -. ^ = 0, (2.6) 

r {rsm9) z 

where the velocities are defined in the form of 3-dimensional contra-variant vectors 

V - dr V - d ° V - d<j) O 7\ 

Vd ~Tv v *~dt- (2 ' 7) 

For the stationary spiral structure globally precessing around z-axis at constant angular 
speed Q (namely the pattern speed), under the coordinate transformation (p = <fi — fit, the 
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solution will be static in coordinate system (t,r,9,<p), namely independent of t. Then we 
get the following equations for such galaxy 

c~ 2 n 2 d 2 $ - [d 2 r + -d r + \(d 2 + cot 9d e + — + AixGp = 0, (2.8) 
(V r d r + Vgd + V v d v )V r - r[V 2 + sin 2 9{V V + fi) 2 ] + <9 r $ = 0, (2.9) 
(V r d r + V e d e + V v d v )V e + *V r V 9 - sin Q cos 9{V V + f2) 2 + ^ = 0, (2.10) 

(V r d r + V e d e + V^)^ + (-V + 2cot9V e )(V tp + fi) + - ^ = 0. (2.11) 

r (rsinf)^ 

In (2.9)-(2.11), some terms have manifest physical meanings, r[V 9 2 + sin 2 0(V^, + fl) 2 } stands 

for centrifugal force, and ( 2 V r + 2 cot 6Vq)VL the Coriolis force. 

Considering an unwarped galaxy with two stationary spiral aims, the total mass-energy 

density and potential can be generally expanded by spherical harmonics Yi m (9, tp) with even 

m. To second order terms, we equivalently have the following approximation 

P = Po + (Pi + P2 cos 2ip + p 3 sin 2ip) sin 2 9, (2. 12) 

$ = $ + ($ 1 + $ 2 cos2<p + $ 3 sin2<p)sin 2 6 l , (2.13) 

where all (p n , $ n ) are functions of r, and (po > 0,pi > 0). Accordingly, the velocity of the 
stars in the disc, also to second order terms, should be 

V r = IVi cos 2^ + W/ 2 sin 2<p, V e = 0, (2.14) 
V^, = co>o + ^i c °s 2(p + u;2 sin 2<p, (2-15) 

where all (W n ,u n ) are functions of r. 

Substituting (2.12) and (2.13) into (2.8), we get the relations between $ n and p n as 
follows 

<9 2 $o + -d r $o + -Ui = ^G Po , (2.16) 

9 2 $! + -a r $i - -^Ui = 47rG Pl , (2.17) 

4 C T 2 fi 2 $ fc + <9 2 $ fc + -d r $ k - = 4vrG Pfc , (A: = 2, 3). (2.18) 

qy ly £ 

Substituting (2.13)-(2.15) into (2.9) and (2.11), and constraining 9 = |7r, we get respectively 
= -d r (Wf + Wi) + <9 r ($ + $i) - ^r(2cj 2 + cj 2 + cj 2 ) + uo x W 2 - u 2 W l + 

2[(U0 Q -Q)W 2 - TUOqUOx + -<9 r $ 2 ] cos 2(p + 

2[{Q - u )Wi - ru u 2 + -d r ® 3 ] sin2<p + • • • , (2.19) 
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and 

= -(rd r u 1 + 2oui)W 1 + -(rd r u 2 + 2uj 2 )W 2 + 

2 

[(rd r u + 2cl> )W / i + 2(cj - Vt)ru 2 + -$3] cos 2ip + 

2 

[(rd r u + 2u Q )W 2 + 2(0 - u )ru! - -$ 2 ] sin2^ + - ■ ■ . (2.20) 

The equation (2.10) automatically holds. Theoretically, we can solve (Q,Wj,Uj) and some 
relations among (/?&, $^) by condition that the coefficients of (2.19) and (2.20) vanish. 

According to the high-accurate observational data[18, 19, 20], we find the rotation curves 
of most spiral galaxies are approximately flat, then we can assume 

u = --n, (r G [Rq, Ri\), (2.21) 
r 

where v is a constant speed with typical value \v\ = 200 ~ 400km/s, [Rq, Ri] is the effective 
region of (2.21), with their typical values as Rq = 100 ~ 500pc, and R\ — 10 ~ 60kpc 
to be approximately the visible radius of the galaxy. Equivalently, we can assume (ujq > 
0, v > 0, Q > 0) in calculation, because the dynamical equations (2.3)-(2.6) have reversal 
invariance under transformation — > —0. We use the empirical condition (2.21) to replace 
the dynamical equations for the background, then the problem is greatly simplified. In 
addition, for a given galaxy, instead of constant v, we can use fitting function v — v(r) in 
(2.21) to get more accurate results and larger effective domain[21, 22, 23]. However in this 
case, we can usually get numerical results only. 

Substituting (2.21) into (2.19) and (2.20), by setting the coefficients of (sin 2ip, cos 2ip) 
terms to zero, we can solve 

ui = -K[lvdr$ 2 + -{v-rn)$ 2 ], (2.22) 
2 r 

oo 2 = -K[\vd r ^ + -(v - rO)$ 3 ], (2.23) 
2 r 

Wi = K[(v-rn)rd r $3 + 2v§ 3 ], (2.24) 
W 2 = K[(v- rtyrd r <$> 2 + 2v® 2 ], (2.25) 
K = (2r 2 tt 2 - ArvVL + v 2 Y l . (2.26) 

By the zeroth order terms in (2.19) and (2.20), we have 

a r ($ + $1) + ^d r (W 2 + W 2 ) = ^r(2u 2 + u{ + cu 2 ) - u x W 2 + w 2 W u (2.27) 
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and 

{rdrUi + 2u l )W l + (rd r u 2 + 2u 2 )W 2 = 0. (2.28) 

Substituting (2.22)-(2.26) into (2.27) and (2.28), we get two constraints for $ fc (r) of the 
background potentials. 

Obviously, the disc satisfies the mass conservation law independent of the dark halo, so 
we have the 2- dimensional continuity equation for all stars and baryonic particles moving in 
the disc as follows [3, 11] 

= <9(£ + V- (V"E), 

= (<9 t + Kdr + + + (2.29) 

where £ stands for the surface mass density of the stars and baryons in the disc. It should 
be mentioned that, (2.8)-(2.11) are expressed in spherical coordinate system due to the 
background gravity, but (2.29) is expressed in polar coordinate system. In the stationary 
case, (2.29) becomes 

(V r d r + VydJE + (d r V r + d v V v + -V r )E = V ■ (VE) = 0. (2.30) 

(2.29) or (2.30) is the equation to describe the mass density of the stars. 



III. RESOLUTION TO THE EQUATIONS 

A. Solution to a Barred Spiral Galaxy 

The general solutions to the above underdetermined equation system are quite compli- 
cated and unnecessary. The most important case to understand the nature of spiral structure 
is the stable and terminal state of a galaxy, which is similar to the eigenstate of a micro 
particle. At first, we consider the case that all stars move in the orbits near circle, for which 
the analytic solutions can be solved. In this case the radial speed is a high order little term, 
then we have 

W x =0, W 2 = 0. (3.1) 

By (2.24) and (2.25), noticing the symmetry between (<3> 2 , $3), we get the equivalent solution, 
except for an initial phase of ip, as follows 

<b 2 = -\(v-rtt)\ $3 = 0, (3.2) 



s 



where q > is a constant. Substituting (3.1) and (3.2) into the above equations (2.16)-(2.18) 
and (2.22)-(2.28), we finally get 

iu x = ^(v-rtt), u 2 = 0, p 3 = 0, (3.3) 
P2 = ((3r 2 fi 2 - 6rvn + 2v 2 ) - ^> - rfi) 2 ) , (3.5) 

2 

= $ + $ 1 + J—(Q r 2 fi 2 - SvrQ + 3v 2 ) - v 2 In — , (3.6) 
24r 4 r 

where r > is a constant with length dimension. 

Substituting the solutions into (2.9) and (2.11), we can check the truncation error of the 

equations reads 

2 2 

A(2.9) = -^r(w-rfi) 2 cos4^, A(2.11) = ~(v - rfi) 2 sin4<p, (3.7) 

which are higher order terms. So the above solution is a good approximation. If q = 0, the 
solutions are exact, which correspond to the galaxy without spirals in the disc. To get more 
accurate results, one should introduce higher order terms such as (cos 4<p, sin 4<p) terms in 
(2.12)-(2.15). 

By (2.21)-(3.2), we get the speed of the stellar flow 

V r = 0, V e = 0, V v = - + -|0 - r ty cos 2cp. (3.8) 

Substituting it into (2.30), we get the continuity equation as d^iV^T,) = 0, which yields the 
mass density distribution of the stars as 

Ur, V )= ■ f B{r l - , (3.9) 
vr z + q[v — riZJ cos lip 

where g(r) is a density function determined by boundary conditions. Apparently, the density 
distribution of the stars display a profile of barred spiral galaxy. 

To determine (p ,Pi), we need another condition related to the dynamics of the black 
halo. In the case of p\ oc po, we can solve 

,2 



n 2 - n - 4 ( 9v 2 fl 2 



Po r~F> 9 



AttG \2(n - 4)(n + 3)r 6 (ra - 2)(n+ l)r 
" 2 -"- 4 ^ (3.10) 



4ttG \^ r 2+n n ( n _ l) r 5 

(3-n)(ra-2) 
n 2 — n — 4 



Pi — 2 — I — ~ Po, (All) 
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in which a and + 1) < n < 3 are constants. (3.10) and (3.11) provide a heuristic 

mass density for the dark halo in barred spiral galaxy within the range [-Rcb -Ri]- 

B. Solution to a Spiral Galaxy 

In an ordinary spiral galaxy, the potentials $2 and $3 should take the following form 

$ 2 = P- I cos(er + ^ ), $ 3 = P" 1 sin(er + ^o), (3.12) 

where £ is a constant, ipo the initial phase of ip and P(r) a function of r. Obviously we can 
set ip = by a translation of <£> similar to the barred spiral case. Substituting (3.12) into 
(2.28), we get a linear equation for P 

P --, F^ P + U +-? F^ + — + ^ P = - 3 - 13 

The equation (3.13) is independent of cos(£r + ip ) and sin(£r + ip ), which reflects (3.12) 
touches the nature of the spiral structure. 

The solution to (3.13) can be expressed by the Hankel function with complex parameters, 
which is much complicated. Here we examine the simplest case of (3.13), the static case with 
Q = 0. The analytic solution to this case can be expressed in a clear form, and it is heuristic 
to reveal the nature of a spiral galaxy. In the general cases, the numerical simulation is more 
convenient and efficient. 

Setting f2 = 0, (3.13) becomes a Bessel-like equation 

P" - -P' + (£ 2 + \)P = 0. (3.14) 

The solution to (3.14) is given by 

P = V^[d J a (|e|r) + C 2 J_ a (|£|r)] , a = ^V7, (3.15) 

where (Ci, C2) are constants, (J a , J_ a ) are Bessel functions with complex parameters, which 
are defined by 

J ^ = ®"£h.WTkTT)&^ { ^ ±ay (3 - 16) 

Substituting (3.15) into (2.18) and (2.22)-(2.26), we can check (p2, P3) and speed all have 
spiral structure like (3.12). Again by (2.30), we find E also has spiral structure. Substituting 
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the results into (2.27), we get a constraint for (<3>o,$i) similar to (3.6). Since the concrete 
expressions are long and complicated, and can be derived by straightforward calculation, we 
do not display them here. 

IV. DISCUSSION AND CONCLUSION 

In the context of general relativity, the whole dynamical equation system for the galactic 
evolution should be the Einstein's field equation 

G^ v = R^ v - -g^R = -kT^, («=^), (4.1) 

combined with the energy-momentum conservation law and the equation of state of the 
gravitating source. For the spinors with interactions, the classical approximation gives[16] 

T" v = ( Ptat + P)U'W + (W -P)g' M/ , (4.2) 

where W is a potential corresponding to the interaction terms, which acts as negative pres- 
sure. According to the energy- momentum conservation law or Bianchi identity = 0, we 
can derive the continuity equation U^T^ = and the equation of motion for the source as 
follows 

Wd^pu* + W) = -(p tot + P)U%, (4.3) 
( Ptot + P)WU% = (<T - Um v )d v {P - W). (4.4) 

For the nonlinear dark spinors, we have W ~ p to t ^ P- In this case, by (4.4) we find 
the stream lines of the spinors are quite different from the geodesies U U U^. U = 0. So unless 
the nature of the dark matter is disclosed, a fully relativistic simulation for the dynamics of 
galaxy is impossible. In (2.1), the effects of (P, W) are merged into one effective mass-energy 
density p, so the treatment is much simplified. 

The above procedure provides a method to research the structure of a galaxy via dy- 
namical approach, which naturally connect the hydrodynamics with empirical data. The 
above solutions provide manifest analytic relations and functions for a spiral galaxy, which 
are much helpful to understand the structure and property of the galaxy. These solution 
verifies that the spiral arms are stationary or even static stellar density wave distribution. 

In the case of barred spiral galaxy with Q 7^ 0, by (3.5), we find = has a root 

r = R < v/Q, which means the bar vanishes at radius r = R ~ R 1 . By (3.3), we get 
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u>i(R)=0, which implies the density wave vanishes near r = R, and then a ring of stars will 
form. The above solution is derived under the assumption (2.21), so the effective domain of 
the solution is also [Rq,Ri], the effective domain of (2.21). These conclusions can be used 
to estimate the global angular speed 

Q « ^ ~ 300km/s/(30kpc) = lOkm/s/kpc ~ 0.002"/year, 
Ri 

which is less than the previously estimated pattern speed 30 ~ 60km/s/kpc[24]. 

By (3.7) and (3.8), we learn all stars in the disc move in the almost circular orbits. This 
conclusion is reasonable and coincident with facts. Through some mechanisms and long- 
term evolution, the baryons were shifted and concentrated into the disc, and their orbits 
of motion became harmonious in a regular galaxy. Otherwise violent collisions will occur 
frequently among the moving stars, and then the galaxy become a hell. 

By (3.3) and (3.5), we find uj\ > and p2 > hold simultaneously. This means where the 
mass density of the background is larger, where the potential is lower and the stellar speed 
is higher. On the other hand, by (3.9) we find that the higher the stellar speed, the lower 
the number density. This means the mass density of the disc should be compensatory for 
that of the dark halo. Namely, the density of dark halo p should take minimum in the spiral 
arms where E take maximum or vice versa. This is a strange but interesting phenomenon, 
which reflects the complicated stream lines of the dark halo. This phenomenon might also 
be one reason why the dark matter and dark energy can be hardly detected in the vicinity 
of the solar system. 

The relations (3.1) and (3.4)-(3.6) as well as (3.12) provide some information for the 
distribution and property of the total mass density and potential. How to combine the 
results with the dynamics of the background is an interesting problem, which might be a 
shortcut to study the properties of the dark matter and dark energy. 
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V. APPENDIX: THE DERIVATION OF THE EQUATIONS (2.1) AND (2.2) 

In this appendix, we derive the dynamics (2.1) and (2.2) of the Newtonian gravitational 
system from the Einstein's relativistic dynamics (4.1)-(4.4) by weak-field and low-speed 
approximation. Some fundamental contents can be found in [17], but here we make more 
systematic and detailed study for galactic dynamics. For convenience, we take c = 1 as 
unit of velocity. Noticing the facts that the collisions among stars rarely occur, and the 
trajectories of the ordinary matter such as electrons and baryons are almost geodesies, so 
for the stars, the following zero-pressure and inviscid energy-momentum tensor holds 

T? = p s U»U v , (5.1) 

in which p s is the comoving mass density of the stars, and is the 4-vector speed of the 
stellar flow. Since the ordinary matter satisfies the mass-energy conservation law indepen- 
dent of the dark halo, we have T^ v v = 0. Expressing it in the form of equations of continuity 
and motion, we get the dynamical equations for the stars 

U^p s + p s U% = 0, U V U% = 0. (5.2) 

The total energy-momentum tensor of the galaxy is still given by (4.2), and satisfies the 
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dynamical equations (4.3) and (4.4). By (4.1) and (4.2), we get 

R = K(p tot +AW-3P), (5.3) 

where R = g^R^ is the scalar curvature. Substituting (5.3) into (4.1), we get 

= -«(Aot + P)U»U V + l -K( Ptot + 2W- P)gT, (5.4) 

where is the average 4-vector speed of all gravitating source. 

In order to make weak-field approximation, we choose the harmonic coordinate system, 
which leads to usual Cartesian coordinate system when making linearization of metric. Then 
we have the de Donder coordinate condition 

F" = g^V% = -J-driy/gfT) = 0, (5.5) 

where g = \ det(g)\. Denote the Minkowski metric by 

V, = rT = diag(l, -1,-1,-1). (5.6) 

For weak-field approximation, we have the linearization for the metric 

g,u = V,u + V> 9^=^ ~ bT, (5-7) 
= rTv vP h a p, h = h^ = rTh^, (5.8) 
g = 1 + h 7 y/g=\ + \h. (5.9) 

In what follows, we directly use = to replace =. By straightforward calculation, we get the 
linearization for other parameters 

1 

2' 

r" = d^-^rTh), (5.11) 

R„v = \d a d a K v - \{i] m d u Y a + Vvad^T"), (5.12) 

BT = ^d a d a h^ v -h,^ ia d a T v + r 1 l/a d a T^) 1 (5.13) 

R = ]^d a d a h - d a T a . (5.14) 
In the harmonic coordinate system, we have 

V = d v {hT - \f v h) = 0, (5.15) 

R,u = \da&*K v , R» u = l -d a d a h»\ (5.16) 

R = \d a d a h, GT = \d a d a (W v - \rTh). (5.17) 



Kp = nVT&cKp + dph av - d u h af3 ), (5.10) 



15 

By (5.15) and (5.17) we find if = at any given time t = to, it will always hold due to 
the Bianchi identity = 0. 

In order to compare with electromagnetism and to understand the physical meaning of 
the parameters, denote 

$ = Ihtt = l -h l \ A= (h tx ,h ty ,h tz ) = -(h tx ,h ty ,h tz ), (5.18) 
H = (h ab ) = (h ab ), ({a, b} e {1, 2, 3}), B = V x A (5.19) 

In the International System of Units, we have the order of magnitude for the components 
of metric 

c 2 |M ~ c|A fc | ~ |$| < 1, (a ^6), (5.20) 

which means \h a b\ <C \Ak\ <C |$| <C 1 if taking c = 1 as unit. 
For the present purpose, we define the stellar speed V by 

V = ^(U\U 2 ,U 3 ), (5.21) 

which is approximately equivalent to the usual definition. For galaxies, we have the following 
order of magnitude 

\V\ ~ 300km/s = l(T 3 c, A ~ kV, h ab ~ k\V\ 2 , (a^b), (5.22) 

in which the coefficient k is also a number of little value. Then according to 

1 = y/g^UW = (1 + 2$ - 2A ■ V + g ab V a V b )^U\ (5.23) 

by omitting 0(V 2 ) terms, the low-speed assumption gives 

U° = 1 -& + A-V. (5.24) 

Substituting (5.21) and (5.24) into (5.2) and omitting the high order terms, we get the 
continuity equation and motion equation for stars 

(dt + V- V) Ps = -p.[V ■ V + (<9 t $ + V ■ A)], (5.25) 
(d t + V- V)V = -V$ + (—d t A + Vd t $) + V x B + V ■ d t H. (5.26) 

In (5.25), we used the de Donder condition T = in the form 

\d t {h xx + h yy + h zz ) = + V ■ A). (5.27) 
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The equation of motion (5.26) has a similar structure to the electrodynamics. From it 
we learn that, $ gives the Newtonian gravitational potential, and A leads to gravimagnetic 
field B. By (5.22), the zeroth order approximations of (5.25) and (5.26) are just (2.29) and 
(2.2) respectively. 

By (5.3) and (5.17), we have 

d a d a h = 2K(p tot + AW -3P). (5.28) 

By (5.28), (5.4) and (5.16), we get the dynamical equations for h^ v 

d a d a h^ = -2K(p tot + P)WW + n( Ptot + 2W- P)rT, (5.29) 
d a d a X ^ = -2n[( Ptot + P)WU V + (W - P)rf v \, (5.30) 

where = - \ifh. If the average speed of the dark halo is also small, omitting 0(U 2 ) 
from (5.29) we get h xx = h yy = h zz = 2*, h ab = 0, (a ^ b) and 

d a d a § = ^d a d a h 00 = -AnGp, (5.31) 
d a d a V = l -d a d a h kk = -4nGp, (5.32) 

where p and p are the effective mass density. Their zeroth order approximation gives 

p = p to t[2(W ) 2 - 1] - 2W + P[2(U ) 2 + l]=p tot -2W + 3P (5.33) 
p = Ptot + 2W-P. (5.34) 

The equation (5.31) is just (2.1). For the dark matter or dark energy with large enough 
negative pressure, by (5.33) we may even have p < 0, which means the Newtonian gravity 
becomes repulsive in this case. So detecting the dynamical behavior of the galactic dark 
halo may be a shortcut to investigate the weird properties of dark matter or energy. 



